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Within the infinite series of ring (or bubble) diagram approximation for the electronic self-energy 
as appropriate for the long-range Coulomb interaction, we calculate the density-dependent T = 
Fermi liquid quasiparticle effective mass renormalization in two dimensional (2D) electron systems 
in the random-phase-approximation (RPA). We find that the quasiparticle effective mass increases 
very strongly with increasing r s , exhibiting an unexpected divergence at a critical r s value (~ 16 
in the ideal 2D system). We discuss the possibility for and the consequences of such an interaction 
induced effective mass divergence as the electron density decreases, and critically comment on the 
likely theoretical scenarios which could lead to such a low density effective mass divergence in 
electron liquids. 
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I. INTRODUCTION 

An interacting homogeneous quantum electron sys- 
tem in a positive jellium background is characterized (at 
T = 0) by the interaction parameter r s , the dimension- 
less average separation between electrons measured in the 
unit of effective Bohr radius, with large (small) r s being 
the strongly (weakly) interacting regime corresponding 
respectively to low (high) densities. The interaction pa- 
rameter r s depends only on the density (n) of the elec- 
tron system with r s oc n~ 1 / 3 (n -1 / 2 ) in 3D (2D) systems. 
The r s -parameter is essentially proportional to the di- 
mensionless ratio of the (interacting) potential energy to 
the (noninteracting) kinetic energy, with the potential 
(kinetic) energy varying as (rj 2 ). (This simple one 
parameter picture of interacting electron systems is valid 
only at T — and for clean free-electron jellium systems 
with little disorder.) 

Alkali metals (e.g. Li, Na, K, Rb, Cs) are examples 
of interacting three dimensional (3D) electron systems 
with reasonably simple band structures so that the ho- 
mogeneous free-electron jellium model may be applica- 
ble. The typical r s -parameter for 3D metals is in the 
range of 3 - 6 (e.g. r s = 3.2 [Li], 4.0 [Na], 5.0 [K], 
5.2 [Rb], 5.6 [Cs]), making metals strongly interacting 
electron systems since r s > 1. Electrons (or holes, as 
in p-GaAs systems) confined by external electric fields 
near interfaces in semiconductor structures (e.g. Si MOS- 
FETs, GaAs quantum wells and hctcrojunctions) are ex- 
amples of 2D electron systems, which are actually quite 
ideal free electron jellium systems if the impurity disorder 
is not too high since the noninteracting energy has free 
electron parabolic form without any band structure com- 
plications. One particular advantage of semiconductor- 
based 2D electron systems is that often (as in Si MOS- 
FETs and GaAs HIGFETs) the carrier density can be 
tuned over one or two orders of magnitude (by changing 
the external gate voltage), thus enabling the tuning of 
the r s -parameter in the same sample over one order of 
magnitude - for example, in high quality Si MOSFETs, 



r s could be tuned from 1 to 12 and in gated n-GaAs 2D 
HIGFET structures, where extremely low carrier den- 
sities can be achieved and r s can be made as high as 
15. In hole-doped 2D p-GaAs structures it is possible 
to achieve extremely high values of r s (~ 20 — 30), but 
there are considerable band structure complications in 
2D p-GaAs hole systems due to the presence of strong 
spin-orbit coupling and the complex multi-band struc- 
ture of the GaAs valence band so that 2D p-GaAs sys- 
tems may not be thought of as simple isotropic free- 
electron-like systems. This capability of tuning the r s 
parameter of a single system over a large range has made 
semiconductor-based 2D electron systems attractive sys- 
tems for studying Coulomb quantum many-body effects, 
and indeed many-body renormalization has been actively 
studied in 2D electron systems over the last thirty years 
both experimentally Ill B l ABIa l a lSIMB and theoret- 
ically [TflllJIlllBLLiM 

In this paper we study theoretically the quasiparticle 
effective mass (m*) renormalization in an interacting 2D 
electron system at T = finding an interesting (and in- 
triguing) divergence in the calculated renormalized effec- 
tive mass at low (high) carrier densities (r s ). Our result 
takes on particular significance since such an effective 
mass divergence has recently been claimed to have been 
experimentally observed in 2D Si MOSFETs with r s val- 
ues close to (around r s ~ 10 — 17) yjSui where we find 
our theoretical divergence. Because of the necessarily 
approximate (and somewhat uncontrolled) nature of our 
theory we cannot be absolutely certain that the observed 
effective mass divergence in 2D semiconductor systems is 
indeed the same phenomenon as what we obtain in our 
theory, but it seems likely that the observed divergence is 
related to our theoretical findings because our calculated 
functional form of m*(r s ) is very similar to the experi- 
mental observation. We emphasize that the self-energy 
approximation (the so-called 'GW' or the single-loop cal- 
culation) used in our theory has been extensively used in 
the literature up tor s ~ 6-7, and our main contribution 
is to extend the calculation to even larger values of r s , 
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mainly motivated by the recent experimental claims we 
mentioned above. 

The existing paradigm (and one of the key cornerstones 
of condensed matter physics) for interacting electron sys- 
tems is Landau's Fermi liquid theory developed fifty years 
ago. The Fermi liquid theory asserts a one-to-one low en- 
ergy correspondence between the noninteracting Fermi 
gas and the interacting Fermi liquid with weakly inter- 
acting quasiparticles (with renormalized single particle 
properties such as effective mass, specific heat, magnetic 
susceptibility, compressibility, etc.) in the interacting 
system behaving similar to free electrons in the nonin- 
teracting system. The key idea is that the Fermi surface, 
i.e. the T — discontinuity in the single-particle momen- 
tum distribution function at Fermi wave- vector k = kp 1 
survives the effect of interaction (albeit with a reduced 
or renormalized discontinuity at Uf) leading to qualita- 
tively similar leading-order low temperature behavior in 
the single particle thermodynamic properties such as spe- 
cific heat, spin susceptibility and compressibility, for in- 
teracting and noninteracting systems. The Fermi liquid 
theory provides a natural explanation for the remark- 
able success of the noninteracting Fermi gas theory in 
describing thermodynamic, magnetic and even transport 
properties of metals, which are in fact strongly interact- 
ing Fermi liquids (r s ~ 3 — 6). The question, however, 
naturally arises whether the Fermi liquid theory remains 
valid as the interaction strength is increased by going 
to higher values of r s . This question has taken on sig- 
nificance in view of recent experiments in 2D electron 
systems at large r s values, where the measured quasi- 
particle effective mass is seen to increase very rapidly 
with increasing r s [E IE IS @ ■ This is the motivation for 
our theoretical investigation of the quasiparticle effective 
mass in 2D electron systems as a function of r s going to 
low carrier densities (i.e. high r s ). 

There are several possible theoretical scenarios for the 
breakdown of the Fermi liquid theory in a clean interact- 
ing electron liquid at low densities (or high r s ). At very 
large r s the system could obviously reduce its potential 
energy (at the cost of higher kinetic energy) by becoming 
an electron crystal (the so-called quantum Wigner crys- 
tal [23). The T = transition from an electron liquid to 
a Wigner solid should be a first order transition in a clean 
system, occurring at r s ~ 35 — 40 for 2D systems as de- 
duced from quantum Monte Carlo calculations [24[ ■ The 
other scenario is the possibility of an interaction driven 
phase transition to some other ordered phase (e.g. charge 
density wave, spin density wave, superconductor, ferro- 
magnet). We do not consider the possibility of such in- 
teraction driven phase transitions in the current work, 
assuming that the system remains a homogeneous elec- 
tron liquid for all values of r s we investigate. Our goal is 
to calculate the 2D zero temperature quasiparticle effec- 
tive mass m*(r s ) as a function of the density parameter 
r s in the interacting electron liquid. 

One problem, which has no rigorous solution, that im- 
mediately arises in this context is that the electron sys- 



tem for r s > 1 is by definition a strongly interacting 
system with a formal failure of the usual Feynman-Dyson 
diagrammatic perturbation theory, and in fact there is no 
known controlled technique for systematically calculating 
exact quasiparticle properties in interacting electron sys- 
tems for r s > 1. In the weak coupling limit, r s <C 1, it 
has been known for a long time that an asymptotically 
exact many-body perturbation theory can be developed 
for quasiparticle properties by expanding in the r s pa- 
rameter. In fact, in the r s -C 1 limit, the expansion 
of the electron self-energy (as well as the polarizability) 
in the infinite series of ring or bubble diagrams (the so- 
called 'RPA', the random- phase-approximation) is exact 
because of the long-range nature of Coulomb interaction. 
Thus in the r s — > limit it is possible to obtain exact 
results for the effective mass renormalization in the in- 
teracting electron system (interacting via the Coulomb 
interaction). Such an exact result is of course useless for 
real interacting electron systems (e.g. metals in 3D, semi- 
conductor layers in 2D) which are never in the weakly 
interacting r s — > limit, but are actually in the strongly 
interacting r s > 1 regime. Therefore, analytical results 
obtained by expanding in r s (assumed to be a small pa- 
rameter) are completely irrelevant for any practical pur- 
pose as was already emphasized by Hcdin |25| and by 
Rice [26| a long time ago, since real 2D and 3D electron 
systems are always in the r s > 1 regime. There is actu- 
ally a deeper (and somewhat subtle) problem with the r s 
expansion of quasiparticle properties in the r s —* limit. 
It has been known for a while [271 that the r s -expansion 
for the interacting electron liquid problem formally fails 
at T = even in the high-density r s — > limit since 
the radius of convergence for such a power series expan- 
sion in r s is, strictly speaking, zero. In fact, a weakly 
interacting (r s <C 1) system of fcrmions does not remain 
a normal Fermi liquid down to zero temperature, but 
makes a transition to a superconducting state (arising 
purely out of the repulsive electron-electron interaction) 
in some rather high orbital angular momentum channel 
with an exponentially small value of T c . In our discussion 
we assume that this so-called Kohn-Luttinger supercon- 
ductivity can be ignored since T c is exponentially small 
for this transition. We uncritically assume the existence 
of a normal Fermi liquid for the purpose of this work. 

We therefore need a robust (but necessarily inexact ex- 
cept in the r s — > limit) many-body approximation for 
calculating the quasiparticle effective mass at arbitrary 
values of r s . Following Hedin and Rice [25|,|26| we use the 
RPA as our central approximation although we improve 
upon RPA in some situations by using various approxi- 
mation schemes as described in details later in this paper. 
We do want to make an important point here about the 
validity (or lack thereof) of RPA in the strongly interact- 
ing r s > 1 regime. It is often stated or assumed that RPA 
is invalid for r s > 1 . This assertion is false as has already 
been argued |23, |2|| |28| a long time ago in the context 
of the applicability of RPA to 3D metals (r s w 3 — 6). 
For arbitrary values of r s , it is more appropriate to think 
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of RPA as a dynamical Hatree-Fock self-energy calcula- 
tion (using the dynamically screened interaction) rather 
than a perturbation expansion in r s (which makes sense 
only in the r s — > high density asymptotic limit). Such 
a dynamical mean- field calculation (called the 'GW' ap- 
proximation (2^, where G is the electron Green function 
and W the dynamically screened Coulomb interaction) 
has a long history in electronic systems, and is exten- 
sively employed in quasiparticle band structure calcula- 
tions (going beyond LDA) in 3D systems (metals, insula- 
tors, and semiconductors) with very impressive empirical 
success. The RPA is a systematic expansion in the dy- 
namically screened interaction (i.e. an infinite resumma- 
tion in terms of the bare unscreened Coulomb interaction 
which has the long-range divergence making an expan- 
sion in bare interaction, even on a formal level, impossi- 
ble), and in general, the validity of RPA for arbitrary r s 
is not rigorously known except that it is rigorously known 
to be exact in the r s — > limit. The reverse is, however, 
not true, i.e., RPA is not necessarily a bad approxima- 
tion for r s > 1 (where one needs to do the calculations 
numerically, of course, as we do in this work instead of 
carrying out an r s expansion) . The applicability of RPA 
for r s > 1 (through exact numerical calculations without 
any r s expansion) must be ascertained empirically, and in 
view of its widespread empirical success in both 3D (met- 
als) and 2D (semiconductor space charge layers) electron 
systems in the r s > 1 regime, we have decided to carry 
out an RPA-based quasiparticle effective mass calcula- 
tion in 2D systems for larger r. values in the regime of 
current interest in 2D systems [E IE IE El ■ We do pro- 
vide some qualitative and quantitative justification for 
the possible validity of RPA self-energy calculations at 
large r s in section llVl of this paper. We argue that our 
RPA calculation of m*(r s ) giving a divergent effective 
mass at a large r s value (~ 16 in the ideal 2D system) is 
qualitatively valid although the precise critical r s value 
where the divergence occurs may not be quantitatively 
accurate within RPA. Our conclusion about the quali- 
tative validity of the RPA prediction of a divergent 2D 
effective mass at large r s is further strengthened by our 
finding that the inclusion of higher order vertex correc- 
tions in the theory (both in the polarizability and in the 
self-energy in a consistent manner) does not change our 
qualitative conclusion at all. 

We believe that the terminology 'RPA' is a misnomer 
for the leading-order dynamical screening self-energy cal- 
culation we are carrying out in this paper although this 
has now become the standard terminology in the liter- 
ature. A more appropriate name for this particular ap- 
proximation scheme (often also called the 'GW approxi- 
mation' in the context of electronic structure calculation) 
is the 'Ring Diagram Approximation' or 'RDA', which 
emphasizes the basic physics of keeping only the infi- 
nite series of ring or bubble diagrams in calculating the 
self-energy. The two most crucial non-perturbative fea- 
tures of the interacting electron physics are preserved (in 
fact, emphasized) in this 'RDA': Screening of the long- 



range interaction to eliminate the Coulomb divergence 
in each order and the collective plasmon excitation of 
the interacting system. Inclusion of these two key non- 
perturbative aspects of the physics of the problem makes 
RPA (i.e. 'RDA') a reasonable approximation for arbi- 
trary values of r s as long as there is no quantum phase 
transition to a new ground state. 

It is well-known [29( that the in r s — > extreme high 
density limit, the RPA calculation for the quasiparticle 
effective mass becomes exact giving the following asymp- 
totic formula (in both 2D and 3D) 
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1 + ar s (b + In r s )+0(r 2 s ), 
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where a, b are constants of order unity in both 2D and 
3D. This formula shows that m*(r s ) has a nontrivial r s 
expansion even in the r s — > limit where the ring dia- 
grams of RPA asymptotically dominate the self-energy. 
We emphasize that m*(r s ) < m, i.e. the mass renor- 
malization is negative for small r s according to the high- 
density expansion. Our goal in this paper is to obtain the 
RPA effective mass numerically without resorting to any 
expansion - the only approximation we utilize is the ex- 
pansion in effective interaction, which we will argue later 
in section HVI is not an expansion in r s but in some other 
effective parameter 8 (a function of r s ) that may actually 
be small even for r s larger than unity. We will argue later 
in this paper that the actual expansion parameter in our 
theory is not the physical density parameter r s , but a 
more subtle parameter S = r s (l + cr s )~ l where c » 1 
making the theory exact in the high-density 5 — > (i.e. 
r s — ► 0), but reasonably accurate for r s > 1 also. 

As mentioned above, the low-density (r s > 5) Fermi 
liquid behavior of 2D electron systems has taken on spe- 
cial significance during the last few years due to the ex- 
tensive experimental activity [E El 01 IE E| m the sub- 
ject generically referred to as the two-dimensional metal- 
insulator-transition (2D MIT), where lowering density in 
clean 2D electron systems causes an apparently sharp 
localization transition in the system. The critical den- 
sity for the 2D MIT seems to be rather close (perhaps 
even the same) as the density for the apparent effective 
mass divergence [E EJ ■ It may appear trivially obvious 



that the conductivity a 



2 r/m*, where n/r/m* are 



the carrier density / transport relaxation time / effective 
mass respectively, will vanish if the effective mass m* di- 
verges (therefore causing an MIT at the critical r s value 
where to* diverges), but this is not necessarily true. In 
particular, in an ideal system (i.e. a translationally in- 
variant disorder-free jellium-background electron liquid) 
the effective mass entering the conductivity formula is 
the bare (band) effective mass (and not the renormalized 
quasiparticle effective mass), and therefore a strong in- 
crease in the renormalized quasiparticle mass to*, as we 
obtain in this work, should not affect transport proper- 
ties in a translationally invariant electron liquid. The 
reason for this is very simple: In a translationally in- 
variant system the total Hamiltonian can be divided into 



4 



two parts, a center of mass Hamiltonian and a relative co- 
ordinate Hamiltonian, with electron-electron interactions 
affecting only the relative coordinate part which does not 
affect long wavelength transport by virtue of Galilean in- 
variance. Real 2D systems, even the very high-quality 
(i.e. low disorder) systems where the 2D MIT phenom- 
ena are typically observed, are not ideal and have con- 
siderable impurity disorder in them. It is, therefore, pos- 
sible that the renormalized quasiparticle effective mass 
appears in the formula for the conductivity of real 2D 
electron systems, with the 2D metal to insulator transi- 
tion at low carrier densities somehow being intrinsically 
connected with the interaction induced low density effec- 
tive mass divergence phenomenon since translational in- 
variance does not apply to real systems. But the issue of 
the relevant effective mass entering the transport theory 
(e.g. the low temperature conductivity) in a disordered 
system is extremely subtle, and it is by no means clear 
that the many-body quasiparticle effective mass m*(r s ) 
could be used in obtaining the conductivity. Much more 
work is needed to clarify this issue. 

Our finding of a divergent quasiparticle effective mass 
in the strongly interacting 2D electron system is very 
similar to the corresponding situation |3Ct l3l| in nor- 
mal He-3 which is known to be an "almost localized" 
Fermi liquid [33 with a divergent effective mass arising 
from the strong inter-particle interaction (which is short- 
ranged in He-3). Similar to the situation in 2D electron 
systems 0, H, @ , He-3 exhibits a strongly enhanced, al- 
most ferromagnetic, susceptibility. We believe that this 
similarity in the quasiparticle properties of normal He-3 
and low-density 2D electron gas is not fortuitous, and 
"almost localization" (i.e. very heavy or divergent effec- 
tive mass) and "almost ferromagnetism" (i.e. strongly 
enhanced susceptibility) may actually be the generic be- 
havior of strongly interacting Fermi liquids, provided the 
interaction is sufficiently strong, independent of whether 
the interaction is long-ranged (as in our case) or short- 
ranged (as in He-3). The presence of disorder in 2D elec- 
tron systems may, however, add considerable complica- 
tions to the 2D semiconductor systems not present in 
He-3. It is interesting to note that the strong-coupling 
regime (with large m* jrti) is attained in the electron sys- 
tem at low densities whereas the strong-coupling regime 
with large m* jm is achieved at high densities in He-3 with 
the difference arising from the long-range (short-range) 
nature of the interaction in electrons (He-3). 

There have of course been many other calculations of 
m*(r s ) in both 3D and 2D electron systems using the 
RPA self-energy approximation. None of these existing 
calculations has however been extended to large r s val- 
ues (i.e. low density), and therefore the actual diver- 
gence of the effective mass m*(r s ) at a critical r s (~ 16.6) 
we find in this work has not been discussed before in 
the literature. Motivated by the 2D MIT physics, how- 
ever, there have been several recent theoretical discus- 
sions 0, 113, H3 of a divergent 2D effective mass at 
low electron-density, but the underlying physics of these 



works 0, 153, are different from the strong interac- 
tion physics discussed in our work. Although we believe 
(and argue in this paper) that RPA remains (at least 
qualitatively) a good approximation even at low carrier 
densities, the question remains whether the divergence of 
m*(r s ) we discover within the RPA self-energy calcula- 
tion is real (i.e. qualitatively valid) or an artifact. The 
related question is whether the Fermi liquid theory ap- 
plies near the effective mass divergence, or equivalently, 
what the nature of this mass divergence behavior is on a 
fundamental level (i.e. is it some kind of a quantum phase 
transition, and if it is, what is the nature of the quan- 
tum phase on the other side of the transition beyond the 
mass divergence?). Unfortunately we are not in a posi- 
tion to definitively answer these questions. But we could 
certainly speculate on the various theoretical possibili- 
ties. First, we believe our RPA calculation to be qual- 
itatively valid, i.e. we believe that m*(r s ) does indeed 
increase very strongly at low densities, possibly diverg- 
ing at a critical density. Second, the divergent effective 
mass obviously signifies some kind of a strong-coupling 
transition whose nature cannot be surmised on the basis 
of our Fermi liquid calculations since the basic premise 
of the one-to-one correspondence with the noninteracting 
electron gas breaks down at the transition point. These 
are all important and interesting questions which should 
be further investigated in future work. In this paper, we 
only report the finding of a low-density quasiparticle ef- 
fective mass divergence in 2D electron liquids within the 
RPA self-energy calculation. 

In section [n] we present the formalism of our theoreti- 
cal approach. In section ITTll we provide our RPA results 
as a function of density. In section llVl we discuss going 
beyond RPA in the process providing theoretical justifi- 
cation for our calculation. We conclude in sectionlVlwith 
a discussion. 



II. FORMALISM 

We describe the basic RPA self-energy formalism in 
III Al followed by the discussion of the finite quasi-two 
dimensional layer width correction in III Bl We discuss 
appropriate inclusion of higher-order vertex corrections 
(beyond RPA) in III CI (more on this in section V) , and 
the on-shell or off-shell calculation of the quasiparticle 
effective mass in III Dl 



A. RPA self energy 

RPA effective mass calculation simply means that the 
quasiparticle self-energy is obtained using RPA dynami- 
cal screening approximation. Fig. ^ shows the Feynman 
diagram for RPA self-energy, from which we obtain (see, 
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FIG. 1: Feynman diagram for self-energy in the RPA calcula- 
tion. Solid lines denote the free electron Green's function Go 
and the wiggly lines the bare Coulomb potential. Note that 
formally the RPA self-energy is Gou, where u is the dynami- 
cally screened interaction in the ring diagram approximation. 

e-g-, 0) for the 2D self-energy (h = 1) 

where v q — 2i:e 2 /q is the bare Coulomb interaction be- 
tween electrons (and q, k appropriate 2D wavevectors), 



Go(k,w) 



1 - "f(&0 ^ 

w - £k + it] lu - £ k - ir] 



(3) 



is the Green's function for free electrons with the nonin- 

,2 

teracting energy dispersion £k = ^ — A*, A* the chemi- 
cal potential, and e(k, to) is the dynamic dielectric func- 
tion. We use 7] to denote an infinitesimal positive num- 
ber, and uf{x) the Fermi function. At zero temperature, 
uf(x) = 1 when x < and otherwise. Within RPA, 
we have 



e(k,w) = l-u,n(k,w) 



(4) 



with n(k, u>) the noninteracting electronic 2D polarizabil- 
ity (i.e. the bare bubble in Fig. HJ: 



n(k,w) = 2 



d 2 q n F (£q) ~ n F (£q+k) 



(27T) 2 Cq " £q+k + W ' 

Previous calculation shows jjj] that in 2D we have 



n(k,o;l 



(5) 




(6) 



Note that in Eq. and Eq. © , the frequency to can be 
any complex number, and the branch cuts of the square 
roots are chosen such that their imaginary part is posi- 
tive. 

Due to the difficulty with the principal value integra- 
tion and singularities of l/e(k, ui) along the real axis in 
Eq. J5J, it is advantageous to follow the standard proce- 
dure and detour the frequency integration from the real 
axis to imaginary axis [33j. After choosing the contour 
of frequency integration as in Fig. |21 we consider the in- 
tegration 



j> ^w(q, v)G (q + k, v + w), 



(7) 



where iu(q, v) is any complex function that is analytic 
in the upper and lower half of the complex plane. On 
the one hand, Q equals the integration of the integrand 
along real and imaginary axes since the integration along 
the curve part of the contour T (Fig. [5} cancels out. On 
the other hand, JJJ) equals the residue due to the pole of 
Green's function 

(1 - n F (£ q+k ))0(w - £ q+k )w(q, £ q +k - u-irj) 

+ ^F(£q+k)0(£q+k - W)w(q,£q+k - <J + if)), (8) 



where 9(x) 
setting 



we have 
S(k,w) 



1 when x > and otherwise. Now by 



u>(k, w) 



e(k,w) 



(9) 



d 2 q 
(27T) 2 

d 2 q 



VqUFiZn+k)) 

dv 



where 



E ex (k,^) = - 



^rcs 



(k,w) 



(2tt)2 
V S res 

d 2 q 

(27T) 2 

d 2 q 



v q w(q, i/)G (q + k, v 



2ni " 
ti line 



(10) 



(11) 



(2n) 2 

(1 - n F (^q+k))0(uJ - £ q +kMq, Cq+k -u-irj) 
+"F(Cq+k)^(Cq+k - w)«>(q, £ q+k - w + irj) , (12) 

dine/,- .a _. f f dv w ^ iv ) 



! (k,c) 



(27r) 2 7 27T iv + u> — £ q+ k 



(13) 



Note that the static exchange part S ex is separated from 
the dynamical self-energy to improve the numerical con- 
vergence of the line part E^ me . The exchange (or Fock) 
self-energy is obviously frequency independent (but de- 
pends strongly on the wave- vector). There is no Hatree- 
contribution to the self-energy in a uniform electron liq- 
uid. 

It is easy to see that the real part of the self energy 
£', which we need in order to calculate the renormalized 
effective mass, can be written as 



E'(k,Cj) = - j -^^VqUF^q+k) 

d 2 q 



(2n) 2 
d 2 q 
(2n) 2 



w<X(q,£ q+k - ^)[«s(Cq+k + riF^q+k)] 

27T V Z + [U> - £q+k) 2 



where ub{x) is the Bose function at zero tempera- 
ture, (i.e., nsix) — — 1 when x < and otherwise) 
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FIG. 2: (Color online) Contour in the frequency v plane. The 
cross denotes the possible position of the pole of Green's func- 
tion. 



w'(q,£ q +k - w) = Re iu(q, £ q +k - w), and u?(q,w) is 
real since w*(q, = w{c[,—iv) = iu(q,«i/). For com- 
pleteness, we also mention that the imaginary part of the 
self-energy comes only from the residue part of the self- 
energy £ res (k, oj). The imaginary part of the self-energy 
does not play any role in the quasiparticle effective mass 
calculation since this is purely a Fermi surface property. 



B. Finite width correction 

The real 2D systems have a finite layer width in the 
transverse direction which reduces the strength of the 
2D Coulomb interaction. In order to test whether our 
result is model dependent, we include finite width cor- 
rection in our calculation to see how the result changes 
compared to our ideal 2D results. The origin of this 
correction comes from the fact that any 2DES confined 
within a quantum well or heterostructure is physically 
a three dimensional system, with the width in the con- 
fined dimension small but finite. Due to the finite width 
of the quasi-2D electron gas system (such as Si MOS- 
FETs and GaAs/AlGaAs heterojunctions), the Coulomb 
interaction should not be written as the ideal 2D form 
V q = 2ire 2 /q, but should have a finite-width quantum 
form factor F such that pj 



nq 



(15) 



where 

F{x) = (1 + ^8 + 9^ + 3^ + 1 



K sc 8(1 + x) 3 



with 



b = 



( K sc^ 



2 \ 1/3 



\ 487rm z e 2 n* 



K sc '2(1 +x) 6 ' 
(16) 



(17) 



denoting the width of the quasi 2D electron g and 
K; ns are the dielectric constants for the space charge layer 



and the insulator layer, R = (K S c + Ki ns )/2, m z is the band 
mass in the direction perpendicular to the quasi 2D layer, 
and n* = with ridepi the depletion layer charge 

density and n the 2D electron density. We choose na e pi to 
be zero in our calculations. When including finite width 
correction in our calculation, we replace v q by Eq. Ijl5|l 
in both Eq. (J2J) and Eq. J3J. Note that we have chosen 
the finite width correction for a triangular potential con- 
finement as in GaAs heterostructures or Si MOSFETs. 
In Si MOSFETs system, 



n sc = 11.5 (Si), 
m = 0.19m„, 



K ins = 3.9 (Si0 2 ), 
m z = 0.916m e , (18) 



while in GaAs heterostructures, 



k sc = 12.9 (GaAs), 
m = 0.07m e , 



k^s = 10.9 (AlGaAs), 
m z = 0.07m e . (19) 



C. Vertex correction 

Even though there has been little systematic work go- 
ing beyond RPA, it is still possible to include in some ap- 
proximate way the effects of vertex corrections to E(k, lu). 
The purpose here is again to test whether our numer- 
ical result is highly model dependent. In some of our 
calculations we therefore approximate vertex corrections 
to the self-energy by including local-field corrections of 
the Hubbard type in the self-energy and in the dielectric 
function dEJUlEi]- In this approximation the self- 
energy S(k, u>) in Eq. (J2J is modified by a vertex term 
7(k,w): 



E(k,w) = - 



d 2 q f dv v q ry(q,uj) 



(2tt) 2 J 2-rti e(q,i/) 
■ G (q + k, v + uj), 



(20) 



with the dielectric function now also including the same 
vertex (local field) correction: 



;(k,u;) = 1 - « g II(k,w)7(fc,cj). 



(21) 



It is important to emphasize that the vertex correction 
must appear both in the self-energy and in the dielectric 
function for the sake of consistency. Within RPA, 7 = 1 
by definition, and in general 7(k, u>) is unknown. We 
choose 



1 



(22) 



where G(k) is a static local-field correction given in the 
2D Hubbard approximation by (2l| 



G(k) 



2^/k 



2 ^k\ 



(23) 



If the inclusion of this local-field vertex correction sub- 
stantially modifies our RPA results, then our qualitative 
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conclusions about the effective mass divergence become 
highly suspect. We will see this not to be the case. Note 
that in section IfV CI we introduce another form of the lo- 
cal field factor G{k) (and refer to this approximation as 
"HA2" , in contrast to the G(k) factor given in Eq. 
which we call "HA1" approximation) in order to compare 
with our RPA results and discuss the qualitative validity 
of our effective mass results. 



D. Effective mass in on-shell and off-shell 
appr oximat ions 

Even after an accurate calculation of £(k, oS), there 
are subtle issues associated with the calculation of the 
quasiparticle effective mass. The single particle energies 
are given by the positions of the poles of the interact- 
ing Green's functions and are determined by the Dyson's 
equation [^j 



written as 



E(k) 



*1 
2m 



E'(k,£(k)- M ), 



(24) 



where £' is the real part of the self-energy. The effective 
mass of the quasiparticles on Fermi surface can be written 
as 



l-£s'(k,o;) 



in 



1 



ra d yv ( 
T"5k \ 



(25) 



k—kp ,u— 



where kp denotes the Fermi momentum. Note that u) is 
measured from the chemical potential /i, which is renor- 
malized by the electron-electron interaction, /i can be ob- 
tained through solving the following self-consistent equa- 
tion 



kp 
2m 



E'(£; f ,lj = 0; //), 



(26) 



where E'(k, w;/Lt) is given by Eq. i|14[l . with chemical 
potential fi an explicit variable with the noninteracting 
chemical potential being /i = Ep = kp/2m, the 2D Fermi 
energy (i.e. £' = 0). Plugging the result of renormalized 
chemical potential /i into Eq. Ij25(l , we can then calculate 
the effective mass. We name this way of effective mass 
calculation the "off-shell approximation" , in contrast to 
the "on-shell approximation" which we describe in the 
following. 

If the self-energy is calculated only in the lowest order, 
it is not sensible to solve Eq. I|24|l but rather take its first 
iteration so that different orders are not mixed in the 
perturbation theory: 



£(k) = — + E'(k,£ k ). 



(27) 



This approximation has been used and justified in the 
literature [H |U |H |H , and here we call it the "on- 
shell approximation". The effective mass can then be 



m 
m 



1 + f J^'(k,60 



(28) 



k—kf 



where we use the un-renormalized (i.e. the noninteract- 
ing) chemical potential Ep in the right hand side of the 
above equation. 

On-shell and off-shell approximations obviously con- 
verge as r s — > 0, where the interaction is weak and 
the higher-order iterations in Eq. I|24|) vanish making 
Eqs. and as well as Eqs. and (gSJ com- 
pletely equivalent. But for r s > 1, where the system is 
strongly interacting, there is a very large difference be- 
tween the two approximations which increases with in- 
creasing r s . The fact that on-shell and off-shell approxi- 
mations give increasingly different (in fact, qualitatively 
different in our case as we will see later in section ITTT^ re- 
sults with increasing r s calls for a careful examination of 
the theoretical issue underlying these two alternate meth- 
ods for calculating the quasiparticle effective mass (see 
section llVl for more details on this topic) to decide which 
one is the correct approximation to use in a particular 
situation (and why). One could take the easy way out 
and insist that, as often is the case in solid state physics, 
empiricism should decide which of these approximations 
is "better" by comparing with experimental data. (By 
this empirical criterion of agreement with experiment, 
the on-shell approximation certainly wins out as the de- 
cisively better approximation for the 2D effective mass 
calculation since the off-shell RPA self-energy approxi- 
mation for the 2D effective mass produces essentially a 
nearly density-independent effective mass with a rather 
small ~ 10 — 20% mass renormalization for r s > 10, 
which is in sharp contrast with the recent experimental 
claims [a 111 HQ-) We feel, however, that the issue of on- 
shell versus off-shell approximation is sufficiently impor- 
tant to merit serious theoretical considerations on its own 
without resorting to empiricism. Arguments supporting 
the on-shell RPA scheme (over the off-shell one) were ac- 
tually provided in rather compelling terms by Dubois |2S| 
and by Ri ce |2al a long time ago (see also more recent 
references [llj and [lJl for the corresponding 2D argu- 
ments). The basic argument, which we will expand upon 
and discuss in more detail in section llVl is that of mixing 
orders in a systematic perturbation expansion. There is 
no question that if one calculates the exact self-energy 
(an obviously impossible task in our problem) then one 
must solve the Dyson's equation to obtain the interact- 
ing quasiparticle Green's function leading automatically 
to the off-shell effective mass formula (i.e. Eqs. (|24fl and 
The problem arises from the leading-order pertur- 
bative (in the dynamically screened interaction) nature 
of the RPA self-energy (formally E ~ uGo in the RPA 
theory, i.e. E is in the leading order of u and Go), which 
implies that all quantities should be expanded only to the 
leading order in the interaction; otherwise one is mixing 
orders in an inconsistent manner in the theory rendering 
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the theory highly suspect. This means that Eq. I|24|) for 
the quasiparticle energy dispersion -E(k) should only be 
solved to the leading order in £ (i.e. keeping only its first 
iteration, as in the on-shell approximation of Eq. J2ZJ), 
and not to all orders as in the off shell approximation. 
This is explicitly seen by iterating Eq. I|24|l in an order 
by order manner leading to 



E(k) = — + £'(k, &) + p (k, 6 + S'(k, &)) 
+E'(k,£ fc + £' (k,&+£'(k ) &))) + •• 



(29) 



Clearly all the terms inside the square bracket are higher 
orders in the interaction u, and should not be kept (as 
long as £' is obtained in the leading order Gqu RPA ap- 
proximation) for the sake of consistency, leaving us with 
the on-shell approximation of Eq. (|27|l . It is important to 
emphasize that the distinction between off-shell and on- 
shell approximations disappears in the weakly interacting 
r s — ► high-density limit since all the terms inside the 
square bracket are higher-order in r s and disappear as 
r s — > 0. We discuss the issue of on-shell and off-shell ap- 
proximations further in section llVl where the justification 
for (or the validity of) our RPA self-energy approxima- 
tion (at large r s ) is taken up. 
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FIG. 3: (Color online) Calculated 2D RPA effective mass 
within on-shell approximation as a function of r s . Inset: the 
denominator of expression 1281 at different r 3 , which equals 
m/m*. It is clear that the effective mass diverges at r s ~ 16.6 
within the on-shell RPA approximation for the ideal 2D sys- 
tem. 
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III. RESULTS 



A. Effective mass within on-shell approximation 



4- 



We first present the results of our effective mass cal- 
culation using Eq. I|28() . We use the dimensionless 
quantity r s to denote the density of the system, r s = 
me 2 /(fi. 2 y / 7m), showing our calculated results as a func- 
tion of r s where n is the electron density. We concen- 
trate on the high r s (i.e. low electron density) behavior 
of the effective mass. We do mention, however, that we 
have calculated m*(r .) f or all values of r s reproducing 
the earlier results [H El E Q E 13 in the literature 
in the r s < 5 regime. In particular, we quantitatively re- 



produce the analytical 2D behavior 



= 1 



2V2 



In 



s/2 



in the r s — > limit within our numerical calculations. 

From Fig. |31 we can see that the 2D effective mass 
has a divergence at low carrier densities (r s ~ 16.6). 
Examining Eq. H28|) carefully, it is not hard to identify 
the origin of this divergence. The dimensionless term 
JrS]'(k, £k)|fc=fc F m the denominator of Eq. ij^ is a neg- 
ative quantity for large r s , and our numerical calculation 
shows that its magnitude increases as r s increases. At 
r s ~ 16.6, the magnitude of this term passes through 
unity, which means that the calculated effective mass di- 
verges. 

We point out that the effective mass divergence we find 
does not in any way automatically imply a concomitant 
failure of the Fermi liquid theory (of course, we do not 
know what happens at r s values above the critical den- 
sity for the mass divergence). In particular, in Fig. 01 we 
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FIG. 4: Calculated inverse Z factor as a function of r a values. 



present our calculated inverse Fermi liquid renormaliza- 
tion factor Z^ 1 = [1 — <9 w £(k, iv)] k=kF w=0 as a function 
of r s . We see that Z remains well-behaved around where 
m*/m diverges (i.e. Z^ 1 does not manifest any diver- 
gence). Since Z denotes the Fermi surface discontinuity 
at hp, Z means that the Fermi liquid behavior is not 
destroyed. 

To test whether our result is highly model dependent, 
we introduce two kinds of corrections as we explained 
in section III Bl and section III CI Fig. shows the cal- 
culated effective mass including finite width corrections 
with RPA, and the Hubbard vertex correction in the ideal 
2D case. We find that finite width correction suppresses 
the effective mass renormalization, while vertex correc- 
tion enhances it at large r s values. The important point 
is that both of these corrections show clear effective mass 
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divergence - for the finite width correction the critical r s 
for the divergence goes up since the bare interaction is 
softened by finite width effect whereas vertex correction 
actually decreases the critical r s since local field correc- 
tion tends to enhance interaction effects. In fact the in- 
verse mass passes zero as r s goes beyond critical values 
in each of these two cases, and the qualitative behav- 
ior of m*(r s ) is very similar to that shown in Fig. [3] 
Note that in our effective mass calculation with finite 
width correction (results shown in Fig. |3J, we have set 
ft = K sc = K ins an d m z ~ m, which is approximately ap- 
plicable for GaAs hetrostructure systems. We also cal- 
culate within the quasi-2D electron gas picture the effec- 
tive mass of Si MOSFETs as a function of carrier density. 
(These results are shown in Fig. If 21 and discussed in sec- 
tion E|) 
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FIG. 5: (Color online) Calculated effective mass within the 
on-shell approximation at different r 3 values. The solid 
(black) curve is the pure RPA result, the dashed (blue) curve 
is the ideal 2D result with Hubbard vertex correction (HA1), 
and the dot dashed (red) curve is the RPA result with finite 
width correction. Note that all three curves show effective 
mass divergence. In the ideal 2D Hubbard approximation re- 
sult, m* diverges at r a ~ 14.2, while in the finite width RPA 
result, m* diverges at r s ~ 27.8. Note that in the finite width 
correction we have set /tins = k S c = R and m z = m (with 
m being the in-plane band mass), and the correction is much 
stronger than in real Si MOSFET quasi-2D system case as 
shown in Fig. 1121 




50 



FIG. 6: (Color online) Calculated 3D RPA effective mass 
within on-shell approximation as a function of r s . Inset: the 
denominator of expression 1281 at different r s , which equals 
m/m*. It is clear that the effective mass diverges as r s ~ 48 
within the on-shell RPA approximation for the 3D system. 



B. Self consistently calculated chemical potential 



The off-shell approximation effective mass calculation 
requires the chemical potential /j, to be calculated self- 
consistently using Eq. I|2fci|) . 
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It is reasonable to ask whether the low density effective 
mass divergence shown in Fig. is a specific character- 
istic of strongly interacting 2D electron systems or is a 
generic property of interacting electron systems. To an- 
swer this important question we have carried out a similar 
calculation for a 3D electron system interacting via 3D 
Coulomb interaction v q = Aire 2 /q 2 , finding a very similar 
effective mass divergence at large r s (see Fig. |SJ). The 
divergence of m*(r s ) at some large critical r s is thus a 
generic feature (similar to He-3) of strongly interacting 
electron systems with system dimensionality (2D or 3D) 
playing no role. 



FIG. 7: Calculated renormalized chemical potential fi using 
Eq. I26H . Ef denotes the Fermi energy for the corresponding 
non-interacting Fermi system. 



In Fig.|7|we present our numerical result for the renor- 
malized chemical potential. The result shows that the 
chemical potential decreases and approaches as r s gets 
larger. We can also see that or the inverse compress- 



ibility K 
or r, increases. 



2*. 
an 



approaches zero as density decreases 
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C. Effective mass within the off-shell 
approximation 

Using the renormalizcd chemical potential, we ob- 
tain the off-shell effective mass results, which are pre- 
sented in Fig. |S] Our off-shell results differ qualitatively 
from our on-shell results, and no mass divergence oc- 
curs within the off-shell RPA calculation at any r s val- 
ues. In fact the effective mass decreases at high r s , which 
is in qualitative disagreement with the experimental re- 
sults 0, El, |j, la, lI 13 ■ It is also strange that the calcu- 
lated m*(r s ) in Fig. [S] (i.e., the off-shell RPA) shows ex- 
tremely weak many-body renormalization, typically less 
than 30% even for r s as large as 10 — 20. 

1.32r 
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FIG. 8: Calculated RPA effective mass within off-shell ap- 
proximation using Eq. I2."ili . 

We now provide some rather simple (essentially hand- 
waving) arguments to better understand the qualitative 
difference between the results of on- and off-shell approx- 
imations. We examine the r s dependence of the par- 
tial derivatives of the self-energy <9 w £'(k, w)|( fc=fcFj w=0 ) 
and (m/fc)9jfeS'(k,a;)|( )fc=fef , i w=0) appearing in Eq. 
Both of these quantities increase in their magnitude 
with r s . As an oversimplified model, we make the ap- 
proximation that 9 w S'(k, w)|( fc _ fcFi w=0 ) = —ar s and 
(m/k)dkE'(k,w)\{k=k F , w=o) = P r s, where a > /? > 
since our numerical results show that the two partial 
derivatives are opposite in sign and close in magnitude. 
Off-shell calculation for effective mass using Eq. (|25|) 
yields the result (l + ar s )/(l + /3r s ), which increases with 
r s at first saturating at large r s without any divergence. 
On the other hand, the on-shell approximation has both 
of these two derivatives appearing in the denominator of 
the effective mass expression (see Eq. (|28|) . which yields 
1/[1 — (a — 0)r s \. (Note that in this oversimplified argu- 
ment we ignore the difference between the renormalizcd 
chemical potential and non-interacting chemical poten- 
tial.) Now we can manifestly see that on- and off- shell 
RPA approximation give the same m*(r s ) in the r s — > 
limit. However, as r s increases the on-shell effective mass 
keeps on increasing, and diverges at r s ~ l/(a — /3) , 



whereas the off-shell effective mass essentially saturates 
for r s > a(> (3), with a rather small (= a/j3 < 30%) 
renormalization. We have already argued that within 
the leading order "Gqu" RPA self-energy approximation, 
the on-shell effective mass calculation is more consistent 
(and hence, better) than the off-shell approximation. We 
will discuss this issue later in section IIVI again. 

For completeness, we also mention that often in the 
literature the effective mass is calculated using Eq. <|25[) 
and the noninteracting chemical potential rather than in- 
teracting (i (i.e. \i = Ep is used in Eq. (|2"5|l ). This is 
somewhat a mixture of our on-shell and off-shell approx- 
imations as it uses noninteracting \i but self-consistent 
quasiparticle energy, and is logically inconsistent. We 
have also calculated effective mass using this method, 
and the results shown in Fig.[5]are similar to the off-shell 
results of Fig. |S] 
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FIG. 9: Calculated 2D RPA effective mass using Eq. 
while keeping fi — Ef- 



IV. VALIDITY OF RPA 

We have discussed somewhat the validity of the RPA 
self-energy approximation in the context of our quasipar- 
ticle effective mass calculation. In this section we con- 
centrate entirely on this important theoretical issue, and 
discuss it at some depth providing several complementary 
arguments supporting (at least) the qualitative validity 
of RPA even at large r s value of primary interest to us 
in this work. 

We start with some general observations about RPA. 
In the context of the quasiparticle self-energy (and ef- 
fective mass) calculations it implies three distinct steps 
of approximations: (i) Calculate the self-energy in a sin- 
gle loop "Gou" approximation (often called "GW" ap- 
proximation in the literature following Hedin's influential 
publication [25| dating back forty years although the first 
such self-energy calculation was originally carried out by 
Quinn and Ferrell [33]), i.e. in the leading order dynam- 
ically screened interaction; (ii) Obtain the dynamically 
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screened interaction by using the infinite series of bare 
bubble diagrams, i.e. approximate the irreducible po- 
larizability by the noninteracting polarizability "GqGq" 
neglecting all vertex and self-energy corrections to the 
polarizability; (iii) Calculate the effective mass using the 
on-shell approximation based on the rationale that in a 
leading order calculation of the self-energy the quasiparti- 
cle energy should only be calculated in the leading-order 
approximation to the Dyson's (integral) equation so as 
not to mix various perturbative orders arbitrarily and 
inconsistently. Of these three steps in the RPA effective 
mass calculation, the first two are essential approxima- 
tions because going beyond the leading order approxima- 
tion in the dynamically screened interaction and beyond 
the bare noninteracting polarizability ("the bare bub- 
ble" ) in dynamical screening are essentially impossible 
in calculating the full frequency and wave-vector depen- 
dent self-energy S(k, to) necessary for obtaining m*(r s ). 
The third approximation, the on-shell approximation de- 
fined in step (iii) above, is, as explained in section IIIII 
a nonessential approximation which we make because we 
believe it to be a better approximation than the full off- 
shell solution of Dyson's equation. 

We first assert that it is incorrect to think of RPA as an 
expansion in the r s parameter as is commonly done. RPA 
is an expansion in the dynamically screened Coulomb in- 
teraction (i.e. the infinite series of ring or bubble dia- 
grams as in Fig.Q, not a density expansion in any sense. 
For example, the RPA calculation for the plasmon dis- 
persion is exact in the q — > limit for all r s , and in real- 
ity, the RPA-calculated 2D plasmon dispersion provides 
an accurate quantitative description for the experimen- 
tal 2D plasmon dispersion [3 HE HH U P to very low 
densities (r s > 10). (Note that the RPA plasmon disper- 
sion is defined entirely by the polarizability through the 
equation e(q, cj) = 0, i.e. 1 — v q H(q, w) — 0. ) What 
is true is that the RPA self-energy approximation (i.e. 
Fig. ^) becomes exact in the r s — > limit, where the 
quasiparticle effective mass can be written as an expan- 
sion in r s , but the precise expansion parameter even in 
this r s — > limit is not r s , but a = e 2 / (hvp) where Vp is 
the Fermi velocity. Writing the expansion parameter a 
out in terms of r s , we see that it is incorrect to think of 
r s as the relevant expansion parameter for the quasipar- 
ticle effective mass calculation even in the well-studied 
r s -> limit since a w r s /12 (3D), r s /4.5 (2D); and 
therefore the true expansion parameter is much less than 
r s , making the radius of convergence of an r s -expansion 
much larger than unity (perhaps qualitatively explaining 
why RPA works so well for 3D metals where r s ss 3 — 6). 
What is undoubtedly true is that for arbitrary r s (> 1), 
RPA is not exact by any means as it is in the unphysical 
r s — > limit. This statement is, however, not equivalent 
to claiming that RPA is necessarily a very bad approxi- 
mation for large r s . The argument for the failure of RPA 
at large r s is based entirely on the wrong premise that 
RPA is a systematic r s -expansion even at large r s . 

It has long been known, since the seminal work of Gell- 



mann and Brueckner fifty years ago [37|], that RPA (i.e. 
the ring or bubble diagrams) provides the asymptotically 
exact ground state energy for an interacting electron gas 
in the high-density r s — > limit (2i| . This is true both in 
2D and 3D interacting electron systems (but not in ID 
system where interaction effects are nonperturbabtively 
singular and one must use entirely different approaches) 
and leads to a ground state energy expression of the form 
e C orr = ^4 (3) + # (3) In r s (A^ + B^r s lnr s ) in 3D (2D) 
electron liquids in the r s — > limit (where A and B are 
known numerical constants), where e CO rr is only the non- 
trivial interaction or the correlation part of the ground 
state energy per particle (i.e. leaving out the noninteract- 
ing kinetic energy and the Fock exchange energy which 
goes as r~ 2 and rj 1 respectively in both 2D and 3D sys- 
tems). It is important to emphasize that the interaction 
energy is manifestly non-analytic around r s — > (due to 
the In r s term) and therefore RPA is not a simple power 
series expansion in inverse-density even in the high den- 
sity r s — > limit. We emphasize here a surprising (and 
apparently little-known) fact [3^: RPA (i.e. the ring 
diagram approximation) again gives an asymptotically 
excellent result for the ground state correlation energy in 
the extreme low-density r s — > oo limit as well. Although 
this was theoretically demonstrated |3^| by Iwata a long 
time ago, this fact seems to be not very well-known in 
the literature. 



We show in llV Al using a very crude and rather simplis- 
tic power counting arguments, that RPA is most likely 
not an expansion in the r s -parameter for arbitrary r s 
values, but it becomes equivalent to an r s expansion as 
r s — > 0. In section llVBI we revisit the on-shell versus the 
off-shell approximation question for the quasiparticle ef- 
fective mass calculation on a more formal level, justifying 
on a diagrammatic basis the superiority of the on-shell 
approximation over the off-shell one showing that the 
on-shell approximation in fact systematically (albeit im- 
plicitly) includes some higher-order diagrams in the effec- 
tive interaction making it therefore better than just the 
simple single-loop "Gou" approximation. Finally, in sec- 
tion llV O we explicitly calculate some appropriate vertex 
corrections to the 2D self-energy by including different 
simple local field corrections to the bare polarizability 
showing that at least this level of an approximate inclu- 
sion of vertex corrections to the theory does not qual- 
itatively change the 2D effective mass results compared 
with the RPA effective mass calculations, providing some 
additional justification for using the RPA self-energy at 
large r s value. None of our arguments establishes the 
RPA effective mass calculation as a theoretically rigor- 
ous approximation for large r s (and surely it is not), but 
we believe that our arguments are persuasive in demon- 
strating that this is not necessarily a qualitatively wrong 
approximation at large r s . 
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A. The effective expansion parameter S in RPA 

To approximately figure out the effective expansion pa- 
rameter in RPA we reconsider the Feynman-Dyson dia- 
grammatic expansion for the electron self-energy and the 
irreducible polarizability in Fig. 1101 by concentrating on 
the higher-order diagrams left out of RPA. From Fig.llUI 
we see that in each higher order (remembering that the 
expansion is in the screened interaction, not the bare 
Coulomb interaction which diverges in the q — > limit) 



G r i d qdu 



with 



we get an additional factor of 8 

the formal series for the polarizability and the self-energy 
being: 



If =>• GqGq + Gq[GqGqu)Gq + Gq(GqGquGqGqu)G<) 

+ ■■■ , 

£ => G u + G u(G G u) + G u(GqG uGoG u) 

+ ■■■, (30) 

where the first terms GqGq (for IT) and Gqu (for E) are 
what we keep in RPA (see Fig. ^) . 

We now carry out a simple heuristic power counting of 
the factor G G w|^ to get: 
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where we have used simple dimensional forms for Green's 
function Go ~ Ep , 2D screened interaction u ~ 
2ire 2 (kp+qTF)~ 1 ; and phase space d 2 q ~ irkp; duo ~ Ep. 
This leads to 



V2r s 



4tt(1 + y/2r s ) ' 



(32) 



Thus the 2D RPA expansion parameter 8 is essentially 
r s as r s — > 0, but is actually of the form r s /[47r(r s + 
l/v2)]r s ) for arbitrary r s , which is certainly not equiva- 
lent to an r s expansion except in the high-density r s «l 
limit. Thus, this admittedly crude power counting di- 
mensional argument leads to the conclusion that RPA 
remains approximately valid even at large r s , where 8 ~ 
1/(471") is indeed a small parameter. This power counting 
argument holds equally well in 3D systems where 



8 = GqGoii 



r s (l/67r)(4/7r)(4/97r) 1 / 3 
1 + r s (4/7r)(4/97r) 1 /3 



(33) 



using the 3D screened Coulomb interaction and 3D phase 
space integral. Thus, in 3D, 8 ~ r s in the r s — > limit, 
but for r s 3> 1, 8 ~ l/(67r) which is small. (We note as 
an aside that our crude arguments fail completely if there 
is a quantum phase transition to a new ground state at 
some finite values of r s since then the form for Go, u etc. 
will change qualitatively.) 



The crude dimensional arguments of this section, like 
all power counting arguments, are too simplistic to be 
taken rigorously. But, these dimensional arguments 
compellingly demonstrate that, contrary to the very 
widespread popular belief, RPA is an expansion in the 
r s -parameter only in the r s — > limit but not for arbi- 
trary values of r s . In fact, for r s > 1, RPA may very 
well be reasonably well-valid (as in commonly found on 
empirical grounds) because the expansion parameter 8 is 
not necessarily large just because r s 3> 1, but is in fact 
quite small (~ 0.1) numerically. We emphasize that the 
dimensional argument given here is for too simplistic to 
be taken rigorously - our purpose here is to demonstrate 
that at large r s the RPA is not equivalent to an r s expan- 
sion. In particular, one serious flaw of our dimensional 
argument is that it completely fails to produce the lnr s 
terms which are known to be present in the r s — > limit, 
but the very presence of these logarithmic terms in the 
high-density limit underscore the fact that RPA is not a 
simple power series expansion in r s , which is the main 
point we are making here. 



B. Approximate higher-order corrections through 
the RPA on-shell approximation 

In this sub-section we show that the on-shell approx- 
imation to the quasiparticle energy, discussed in sec- 
tion [n] actually implicitly incorporates some higher- 
order corrections (i.e. beyond the single-loop Goit self- 
energy of RPA), making it therefore an improvement 
over the off-shell approximation involving the full solu- 
tion of the Dyson's (integral) equation. The basic ideas 
underlying this argument go back a very long time to 
Dubois |2g and to Rice .26). We outline the arguments 
here in the context of the 2D RPA self-energy calcula- 
tion since the arguments seem not to be very well-known 
(or well-appreciated) in the literature. The iterative (off- 
shell) solution of the Dyson's equation to the quasiparti- 
cle energy can be formally written as fFig. HOf d) ): 



Eu 



£ k + £'(k,£; k ) 

6 t + £ , (k > &) + E'(k,£ k + E'(k,60) 

<9£ 



£ k + £'(k,&) + £'(k,£k) 



(34) 



where we have only explicitly considered up to the 
second-order in iteration. One can now formally com- 
bine the second order iteration in the self-energy in the 
Dyson's equation with the second-order terms in the dia- 
grammatic perturbative expansion (Tig. llOf a)) of E itself 
(i.e. the two second order diagrams in Fig. llOf a) beyond 
the Go« single-loop leading-order RPA self-energy term) 
to obtain the second order correction to the basic 'on- 
shell' RPA self-energy term: 
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where E^ and EI, are respectively the contributions 
from the second order exchange self-energy and the sec- 
ond order vertex correction diagram, and E^ is the first- 
order RPA self-energy obtained from the second-order 
Taylor expansion of the Dyson's equation expansion. 
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FIG. 10: Feynman diagrams for self-energy including higher 
order terms in the perturbation theory which are excluded in 
Fig.E] 

It is now easy to see that the second order exchange 

(2) 

term (T,\ ) and the term arising from the iteration of the 
leading-order self-energy essentially cancel each other: 



,(2) 



0E« 



E« 



d 2 qd 2 q' f dvdv' 



1 d£k J (2tt)4 J (2tt)2 

xG (k + q, £k + is) 2 G (k + q', a + v') 
Xu(q, v) [u(q - q', v' - v) - u(q', v')] . (36) 

By generalizing these arguments order by order, it can 
be shown in a straightforward fashion that such approx- 
imate cancellations between the higher-order exchange 
corrections and the iterations of the Dyson's equation 
for the leading-order (i.e. RPA) self-energy occur to all 
orders. This leads to the conclusion, already made in 
section ITT1 and Hill that the on-shell RPA self-energy is a 
better approximation than the off-shell one since the on- 
shell approximation to the leading-order RPA self-energy 
implicitly contains some aspects of the higher-order cor- 
rections through the order-by-order cancellation demon- 
strated above. In particular, the on-shell 'Gqu' approxi- 
mation simulates the self-consistent 'Gm' approximation 
, where G is the full interacting Green's function rather 
than just the noninteracting Go- 

This reinforces our earlier assertion that the RPA on- 
shell effective mass calculation is not a simple high- 
density expansion, and could have substantial qualita- 
tive (and perhaps even quantitative) validity at larger r s 
values although the extent of its quantitative validity is 
difficult to assess because it is not a controlled perturba- 
tive theory at large r s . 

Similar to the other arguments made in this section, we 
emphasize that the cancellation between the higher order 



diagrams (Fig. Illfl and the Taylor expansion discussed 
here is certainly partial and approximate, and again this 
is by no means a rigorous argument. 



Including vertex corrections through Hubbard 
local field effects 



To obtain some ideas about the model dependence of 
our calculated quasiparticle effective mass we have al- 
ready shown results (see Fig. comparing the 2D RPA 
effective mass with the 2D effective mass calculated in 
the Hubbard approximation. The basic idea is to in- 
troduce an approximate vertex correction (to both the 
polarizability function and the self-energy) in the theory 
to calculate the renormalized effective mass in order to 
assess the robustness of the RPA approximation. A rigor- 
ous calculation of the self-energy vertex correction in the 
dynamical screening expansion is essentially an impossi- 
bly formidable task; one must therefore resort to drastic 
approximations. One reasonably successful technique for 
introducing vertex corrections to the theory has been to 
preserve the formal structure of the RPA self-energy and 
the polarizability with a new function 7, the so-called 
local field correction, modifying the electron-electron in- 
teraction mimicking the higher-order vertex corrections. 
Formally, the self-energy is given by Eqs. I|2(J|I and Q2ip. 
and the 7 function is given in terms of a local field factor 
G(k) by Eqs. and G(fc) can be chosen satisfy- 

ing various sum rules. The two most popular forms for 
the local field correction are due to Hubbard |39j , and in 
2D system they take the form 



HA1 : G(fc) 
HA2 : G(fc) 



k 



2^k 2 + k 2 F + qjp 



(37) 
(38) 



We have already used the first model (HAl) in obtaining 
results in Fig. [5] We find a divergent m*(r s ) for r s ~ 
10 — 17 in all three approximations with HAl giving a 
stronger divergence than the other two. In Fig. ^] we 
compare these two vertex correction approximations with 
RPA for the calculated 2D m*(r s ) as a function of r s 
in the r s = — 17 regime. The results show that (i) 
qualitatively all three approximations give similar results, 
and (ii) the HA2, which is considered 26] to be a superior 
vertex correction than HAl, actually produces m*(r s ) 
quantitatively very close to RPA. 

The results of Fig. 1111 which include vertex corrections, 
again confirm the essential qualitative validity of our cal- 
culated m*(r s ). The point we are making here is that our 
main qualitative conclusion that m*(r s ) diverges at some 
value of r s remains valid even when Hubbard approxima- 
tion the vertex correction is made by going beyond RPA. 
On the other hand the precise value of r s where m* di- 
verges may depend on the approximation. 
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FIG. 11: (Color online) Calculated 2D effective mass within 
the on-shell approximation at different r s values. The solid 
(black) curve is the pure RPA result, and dashed (blue) curve 
is the HA1 result, and the dot dashed (red) curve is the HA2 
result. RPA, HA1 and HA2 results all show effective mass 
divergence. It is clear that the RPA result and the HA2 result 
are strikingly close to each other. 



It is "well-known that the strongly interacting Coulomb 
system at large r s , where the electron-electron interac- 
tion cannot be considered a weak perturbation to the 
noninteracting kinetic energy, does not allow for a con- 
trolled systematic many-body calculation since no order- 
by-order perturbative diagrammatic expansion is possi- 
ble because there is no small parameter (in contrast to the 
high-density situation, r s « 1, where the interaction is a 
weak perturbation and the RPA self-energy expansion is 
asymptotically exact). The great merit of RPA lies in the 
fact that it is a systematic self-consistent mean-field the- 
ory (apart from also being an exact low-r s perturbative 
theory) because it includes the basic physics of the dy- 
namical screening of the long-range Coulomb interaction. 
As such RPA should remain qualitatively valid at large 
r s until there is a quantum phase transition to a new 
state. There have been many systematic attempts |4Cj 
to go beyond RPA in the literature in calculating the 
electron self-energy, and many of these attempts include 
the approximate vertex corrections by exploiting the self- 
consistent field (rather than the perturbative) nature of 
RPA by including the vertex corrections through various 
local field factors. The two Hubbard approximations we 
include in our work are only among the simplest such ap- 
proximations incorporating vertex corrections- we refer 
to the literature 40] for more sophisticated treatments of 
vertex corrections. The important point in the context of 
our theory is that the incorporation of vertex corrections 
is unlikely to cause a qualitative change in the RPA re- 
sults although quantitative changes are expected. There- 
fore, we can conclude that our finding of an effective mass 
divergence within the RPA self-energy calculation (and 
the Hubbard approximation) should remain valid in gen- 



eral although the critical r s where the mass divergence 
occurs is likely to be affected by the specific many-body 
approximation. 



V. DISCUSSION AND CONCLUSION 

We have presented in this paper a detailed quantitative 
calculation for the zero-temperature quasiparticle effec- 
tive mass in a strongly interacting homogeneous electron 
liquid using the realistic long-range Coulomb interaction 
and the noninteracting parabolic kinetic energy disper- 
sion. The theory covers a large range of the dimensionless 
interaction parameter space, and provides m*(r s ) down 
to the low-density regime r s ~ 15 — 20. We have focused 
primarily on the (semiconductor-based) 2D electron sys- 
tems since rather high (low) values of r s (density) are now 
routinely achieved experimentally in 2D semiconductor 
systems, and therefore there is considerable current in- 
terest in the many-body properties of 2DES. Our main 
results for m*(r s ) are qualitatively similar in 2D and 3D 
systems. We also provide in this paper a serious and sub- 
stantive theoretical analysis of the feasible many-body 
approximations for strongly interacting electron liquids 
within the Feynman-Dyson diagrammatic perturbation 
theory. We argue (based on a number of approximate 
non-rigorous arguments) that the RPA self-energy ex- 
pansion in the dynamically screened Coulomb interac- 
tion, the dynamically screened single-loop Hatree-Fock 
(or equivalently 'Gow' or 'GW) approximation, may very 
well be a reasonably good (albeit somewhat uncontrolled 
in the sense that a systematic improvement on the single- 
loop 'Gou' approximation is difficult) approximation even 
at large r s (while being exact in the r s — > limit, as has 
been known for a long time). In the process we estab- 
lish that the RPA self-energy or effective mass calcula- 
tion is not an r s expansion at arbitrary r s values, it is 
more likely to be an expansion in an effective parameter 
8 = (r s /c)/(a + r s ), where a ~ 1 and c > 1 in both 2D 
and 3D. Thus, RPA is exact in the r s — > limit, and is 
reasonably well- valid for r s 3> 1 provided c> 1; we find 
c w 4tt (2D), 6?r (3D) showing that the RPA self-energy 
and effective mass calculation may have reasonable qual- 
itative (perhaps even quantitative) validity in the large 
r s regime of our interest (provided that the strong in- 
teraction does not drive the system to a quantum phase 
transition) . The precise level of the quantitative validity 
of RPA cannot, however, be ascertained at large r s values 
since we have no reliable way of estimating the terms be- 
ing left out of the calculation. (In this respect RPA shares 
the shortcomings of other well-known electronic theory 
approximations such as LDA or DMFT which are very 
successful quantitative theories without being controlled 
systematic approximations.) Analyzing in some depth 
the structure of the infinite resummation of perturbation 
terms involved in the RPA expansion in the screened in- 
teraction for the self-energy function, we also establish 
rather compellingly that the on-shell RPA quasiparticle 
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energy (and effective mass) calculation is the appropri- 
ate correct approximation at arbitrary r s values, and the 
off-shell approximation in fact leads to qualitatively in- 
correct results (by virtue of the inconsistent mixing of dif- 
ferent orders) at large r s ; for small r s (<C 1), of course, the 
on-shell and the off-shell approximations become equiv- 
alent. 

We emphasize that, in spite of several crude (and at 
best persuasive, by no means rigorous and compelling) 
arguments we have provided in section IIVI in support 
of our approximation scheme, our results can only be 
considered qualitative since the ring-diagram approxima- 
tion inherent in RPA is exactly only in the high density 
(r s — > 0) limit whereas our finding of a low-density ef- 
fective mass divergence occurs in the strongly interacting 
r s w 10—15 regime. While agreeing that RPA is, in prin- 
ciple, an uncontrolled many-body approximation at large 
r s , we claim that there is no particular value of r s at (or 
above) which the RPA self-energy fails qualitatively. In 
fact, the evidence is quite contrary, typically RPA predic- 
tions for qualitative behavior always turn out to be valid. 
For example, RPA energetics for interacting electron sys- 
tems predict a ferromagnetic instability at r s w 19(7) 
in 3D (2D) electron systems where the spin <?-factor (or 
equivalently, the static susceptibility) diverges, and it is 
now universally accepted that an interacting electron sys- 
tem indeed undergoes an interaction-driven low-density 
ferromagnetic instability ( the so-called Bloch ferromag- 
netism) albeit at r s values larger than the RPA predic- 
tions. We see no reason why our finding of a low-density 
effective mass divergence within RPA should be any dif- 
ferent from the corresponding divergence in the magnetic 
susceptibility. It is likely that RPA underestimates the 
critical r s where the effective mass diverges, and the real 
effective mass divergence would occur at higher r s values, 
but we see no reason for RPA to give a qualitatively incor- 
rect answer for the effective mass divergence. An impor- 
tant related point in this context is that the RPA or the 
"GW approximation" is the only available many-body 
approximation for quantitative electronic self-energy cal- 
culations, and as such it is routinely used for 3D metallic 
systems with r s «5. In 2D electron systems, RPA effec- 
tive mass calculation has earlier been carried out in the 
literature up to r s m 5 also in the context of providing a 
quantitative explanation for effective mass measurements 
in 2D semiconductor structure. We have now extended 
these calculations up to r s w 10 — 15 because recent low 
density 2D systems have achieved such large values by 
providing high-quality low-density 2D structures. 

Our theory involves two distinct approximations: RPA 
or ring-diagram approximation for the electronic self- 
energy calculation and the on-shell approximation for 
calculating the quasiparticle energy leading to the quasi- 
particle effective mass. While the RPA self-energy cal- 
culation is an essential (and in some sense, uncontrolled) 
approximation in our theory (except for incorporating 
vertex corrections through crude local field approxima- 
tions, e.g. the Hubbard approximation, as we do in some 



of our calculations) , one could question the need (or even 
the validity) of the on-shell quasiparticle effective calcu- 
lation. We have argued in this paper that for our 'Gqu' 
self-energy approximation (Fig. QJ, the on-shell effective 
mass calculation is, in fact, the appropriate approxima- 
tion (rather than the off-shell approximation which solves 
the full Dyson's equation iteratively) . We assert that, as 
shown by Rice [26j and used by Quinn and collabora- 
tors 01 a l° n g time ago, the effective mass approxima- 
tion consistent with the Landau Fermi liquid theory and 
our 'Gou' self-energy approximation is precisely the on- 
shell approximation. This is simply because we use Go, 
the non-interacting Green's function in all our internal 
propagator, and therefore the quasiparticle energy con- 
sistent with this leading-order self-energy calculation is 
precisely our Eq. 1)27(1 . the on-shell approximation, rather 
than our Eq. (|24|l . the off-shell approximation. This has 
been discussed in great details by Rice j2|| almost forty 
years a go ( and was already implemented for 2D systems 
in Ref. [13 almost thirty years ago). Our on-shell quasi- 
particle effective mass is therefore precisely the Fermi 
liquid quasiparticle effective mass 411. 

We find that the quasiparticle effective mass m*(r s ) 
at T = increases rapidly (both in 2D and 3D) at low 
densities with increasing r s , consistent with the strongly 
interacting nature of the Coulomb system at large val- 
ues of r s . One most important and interesting result is 
a true divergence of m*(r s ) at a critical r s ~ 16.6 in 
2D (the mass divergence occurs around r s ~ 48 in 3D) 
within the on-shell single-loop Gqu RPA self-energy cal- 
culation. This theoretical effective mass divergence (or at 
least the rapid increase of effective mass with decreasing 
density) seems to superficially simulate recent measure- 
ments [a, IS 13 °f low density 2D effective mass (from the 
temperature dependence of the SdH oscillations of the 
low field magnetoresistance) in n-Si MOSFETs. 

There are three general questions that immediately 
arise in the context of our theoretical finding of a di- 
vergent m*(r s ) for sufficiently large (but experimentally 
accessible in 2D, but not in 3D systems) values of the 
density parameter r s . These questions are: (i) Are our 
theoretical results of any relevance to experiments (par- 
ticularly to the large body of 2D M-I-T literature where 
large- r s behavior of 2DES is being explored extensively)? 
(ii) Is our finding of a divergent effective mass at a large 
critical value of r s real? (i.e. is it a truly robust find- 
ing as we have argued in this work or just an artifact 
of our single-loop Gqu approximation scheme?) (iii) If 
the interaction-driven divergence of m*(r s ) is real, then 
what does this phenomenon imply for an electron liquid? 
Below we discuss some partial answers to (and educated 
speculations on) these questions. Since an exact solution 
to this strongly interacting (and 50-year old) many-body 
problem is unlikely to come in the near future, we can 
only discuss various aspects of this problem qualitatively 
with no definitive resolution of these important questions. 

For whatever it is worth we can, of course, carry 
out a direct comparison with the experimentally mea- 
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sured density-dependent effective mass with our calcu- 
lated m*(r s ) = m*(n) for the appropriate 2DES. We 
show our theoretical results for the T = RPA effective 
mass as a function of the carrier density in Si-Si02 inver- 
sion layer in Fig. 1121 including the quasi-2D form-factor 
effects in the Coulomb interaction by modifying the bare 
Coulomb interaction v q — 2ire 2 /(Rq) to v(q)F(qb) where 
R is the effective background dielectric constant for the 
system, 'F' is the quasi-2D form-factor, and b = b(n) 
is the finite quasi-2D layer width parameter discussed in 
section III Bl For Fig. ^| results we have used the ap- 
propriate semiconductor electron band mass and lattice 
dielectric constant (given in Eqs. I|18fl ) in the calcula- 
tion. In Fig. 1121 we have shown as an inset the experi- 
mental Si MOS effective mass values for a comparison 
with our theory. Clearly, the theory and experiment 
are in excellent qualitative agreement, most notably in 
the rapidly rising low density m*(n) behavior. While 
it is certainly gratifying to see that our theory agrees 
well with the existing low-density effective mass mea- 
surements, one must take such agreements with more 
than a grain of salt. The reliability and/or the level of 
accuracy of the experimentally measured m*(n), partic- 
ularly at the low carrier densities of interest to us, is 
unknown. This is particularly true in view of the ex- 
perimental effective mass measurements being based on 
the temperature-dependent amplitude of the weak-field 
SdH magnetoresistance oscillations at each density. We 
have shown elsewhere that the quasiparticle effective 
mass m*(r s ) has a strong temperature dependence at low 
densities - in fact m*(T/T F ) ~ m*(0)[l + c(T/T F )\ with 
c > 0. Such an intrinsic temperature dependence of the 
effective mass itself may introduce unknown errors in the 
experimental mass measurements. In addition, there is 
always the key problem of not quite knowing what the 
electron temperature really is (a particularly severe prob- 
lem at low carrier densities) which could lead to larger 
errors in the experimental determination of m*(n). In 
spite of these caveats the excellent qualitative agreement 
between theory and experiment demonstrated in Fig. 1121 
is quite impressive. 

In addressing the issue of the "robustness" of our find- 
ing of a divergence of m*{r s ) at a large critical r s we 
first assert that within our approximation scheme this 
interesting result is certainly real and not an artifact as 
demonstrated by the manifest model independence of our 
result - we find that m*(r s ) diverges at a large r s both 
in 2D and 3D systems, both within the simple RPA and 
the Hubbard approximation calculations, and both in the 
strictly 2D and quasi-2D systems. The only model de- 
pendence of our theory is in our finding that the off-shell 
approximation (i.e. the full solution of the Dyson's inte- 
gral equation using, however, only the leading order dy- 
namically screened single-loop self-energy function) does 
not leads to a mass divergence (and in fact lead instead 
to an effective mass saturation at large r s in clear dis- 
agreement with experimental results) , but we have given 
convincing arguments that the off-shell approximation is 
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FIG. 12: Calculated RPA effective mass within the on-shell 
approximation at different density values for 2D Si MOSFET 
system. Inset: Experimental effective mass results with data 
taken from Ref. 



theoretically unjustified ( and inconsistent) in view of the 
single-loop leading-order nature of our 'Goit' self-energy 
calculation. 

Although the divergence of m*(r a ) at large r s is cer- 
tainly a real result within our 'Gow' single-loop self- 
energy approximation scheme, the important question 
of whether it is a correct result (i.e. valid beyond our 
approximation scheme) for the interacting electron liq- 
uid remains unanswered at this stage. While we have 
included some approximate vertex corrections (see sec- 
tion IIV|) in the theory, a self-energy calculation going 
beyond the leading-order 1 Gqu' single-loop diagram is 
essentially impossible. The second-order self-energy di- 
agrams in the dynamically screened Coulomb interac- 
tion expansion (Fig. I1U|) involve six-dimensional highly- 
singular (with overlapping singularities and branch cuts) 
integrals in the 'GouGoGoii' second-order terms. We see 
no hope of an accurate (numerical) evaluation of these 6- 
dimensional highly singular principal value integrals in 
order to reliably calculate S(k,w) to the 2-loop level 
so that m*(r s ) can be calculated correctly (formally) 
to second order in the dynamically screened interaction. 
Such a (hopeless) calculation is necessary to see whether 
the divergence of m*(r s ) is truly a feature of the in- 
teracting electron systems at low carrier densities or is 
merely an artifact of the single-loop self-energy formal- 
ism. Even such a 2-loop calculation will not obviously 
definitively establish the existence of a density-driven ef- 
fective mass divergence transition (because one could still 
wonder what would happen at the 3-loop or the 4-loop 
level), but we believe that finding a divergence of m*(r s ) 
in both 1-loop (as we do here) and 2-loop level will go a 
long way in establishing the very strong possibility that 
such a divergence does indeed occur in the ground state 
of an interacting quantum Coulomb Fermi system. 

Short of such a 2-loop self-energy calculation (which 
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we believe to be impossibly difficult) we have to make 
do with a number of heuristic arguments supporting the 
reality of the effective mass divergence reported in this 
work. First, the single-loop G$u RPA self-energy cal- 
culation (along with its simple Hubbard approximation 
vertex correction) certainly leads to an effective mass di- 
vergence. Second, we give several arguments (see sec- 
tion IIV|) why we believe the Gqu approximation to be 
at least qualitatively valid at large r s where the effective 
mass divergence occurs. Third, the structure of the ef- 
fective mass formula (see section [DJ indicates that the 
necessary and the sufficient condition for the divergence 
of m*(r s ) is to have ^ — — 1, which is certainly allowed 
as a matter of principle. Fourth, the fact that we get 
the effective mass divergence in both 2D and 3D sys- 
tems indicates that the effective mass divergence may be 
a generic feature of interacting electron systems in the 
strong interaction regime. Fifth, as mentioned in the In- 
troduction, the strongly interacting normal He-3 (with 
strong short-range interactions) is also thought to be an 
almost-localized Fermi liquid (i.e. very large quasipar- 
ticle effective mass) akin to what we find for strongly 
interacting electron liquids interacting via the Coulomb 
interaction. Similar to our finding in this work, normal 
He-3 is also reported [23,ISE3| to go through an effective 
mass divergence transition at sufficiently large interac- 
tion strength. It is therefore quite possible, perhaps even 
very likely, that both short-range interacting He-3 and 
long-range interacting electron liquids have quasiparticle 
effective mass divergence at sufficiently large interaction 
strength; such an effective mass divergence leading to "an 
almost-localized Fermi liquid" may very well be a generic 
feature of strongly interacting quantum Fermi systems. 

Accepting that there is indeed a critical (large) value of 
r s = r c (ss 16.6 in the ideal 2DES according to the single- 
loop calculation, but in general we expect r c to be some 
large r s value indicating the strong interaction regime) 
where m*(r s )\ r ^ r diverges, the important question of 
the nature and implications of this mass divergence takes 
on great significance (both theoretical and practical since 
2DES with r s > 10 is routinely available with the state 
of the arts semiconductor technology) . We first note that 
according to the Landau's Fermi liquid theory the low- 
temperature thermodynamic properties such as specific 
heat, spin susceptibility, and compressibility are all pro- 
portional to the quasiparticle effective mass: 
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where Aq, Bo, A\ are the appropriate Landau Fermi 
liquid parameters. Assuming that Aq and Bq are well- 
behaved at low densities, we conclude that both the spin 
susceptibility and the compressibility will also diverge 
around the critical r s where the effective mass diverges. 



Such anomalous behavior in the magnetic susceptibility, 
with rapidly increasing susceptibility at low densities, 
has indeed been reported in several recent experimen- 
tal studies 0, 13 of 2DES. Recent 2D experiments also 
report |42] | the inverse compressibility passing through 
zero at low carrier densities indicating an infinite com- 
pressibility consistent with the mass divergence (assum- 
ing regular density dependent behavior of the Landau 
parameter 1 + Aq). But the compressibility measure- 
ments are most likely connected with disorder effects 
invariably present in the system, and are unlikely to be 
intrinsic Fermi liquid renormalization phenomena. More 
work is needed to settle these issues. 

A key question in this context (accepting the mass di- 
vergence to be real) is whether this effective mass diver- 
gence signifies a new strong coupling fixed point of the 
interacting electron system or is a subtle manifestation 
of one of the already known quantum phase transitions 
which may occur in an interacting electron system. A re- 
lated question of equal importance is whether this mass 
divergence necessarily implies a breakdown of the Fermi 
liquid theory. Unfortunately, none of these questions can 
be answered at our current stage of understanding of the 
subject since our calculation only leads to the effective 
mass divergence without telling us the nature of the new 
ground state, if indeed there is a quantum phase transi- 
tion associated with the mass divergence phenomenon. 
Using understanding |32| developed in the context of 
the strongly interacting neutral He-3, which also under- 
goes |3(J,|3ll a very similar effective mass divergence tran- 
sition, we speculate that the transition associated with 
the effective mass divergence may very well be a con- 
tinuum version of a Mott transition, where increasing 
Coulomb interaction in a system leads to "localization" 
and metal-insulator transition as the system can lower 
its interaction energy at the "cost" of kinetic energy by 
developing infinite mass (i.e. "no hopping" ) - a divergent 
mass in a continuum system is the closest analog to the 
band- width collapse due to the disappearance of hopping 
in the Mott transition in a lattice system. It is gener- 
ally assumed that the Wigner crystallization transition, 
which in 2D (3D) systems is thought |24Ll43j to occur at 
large r s values ~ 38(72), is the continuum analog of the 
Mott transition, but it is, in principle, possible for the 
system to undergo a different type of localization transi- 
tion through the effective mass divergence long before the 
Wigner crystallization low-density point is reached. The 
alternate possibility is that our effective mass divergence 
is essentially a signature or a precursor to the eventual 
Wigner solidification transition in the system. We also 
can not rule out the possibility that the RPA effective 
mass divergence is a precursor of a charge density wave 
transition in the system (related to, but not identical to, 
the Wigner transition). Any of these transitions of course 
signifies a breakdown of the normal Fermi liquid picture 
because the one-to-one correspondence with the nonin- 
teracting electron gas ground state is destroyed beyond 
the transition point. At this stage, however, we do not 
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rigorously know the nature of the interaction-driven ef- 
fective mass divergence phenomenon (except that it hap- 
pens at some large r s values). Much more work will be 
needed to develop an understanding of the mass diver- 
gence phenomenon reported in this work, but the very 
real exciting possibility exists that all strongly interacting 
homogeneous Fermi liquids undergo a continuum version 
of a Mott transition, where the renormalized quasiparti- 
cle effective mass diverges, at sufficiently high interaction 
strengths. At this stage we do not know whether our 
theoretical discovery of a low-density divergence in the 
quasiparticle effective mass liquefies any kind of phase 
transition or not. In may be worthwhile to point out 
in this context that the corresponding problem of a sin- 
gle electron strong coupled to the lattice ("the strong 
coupling polaron" problem) also shows a sharp effective 
polaronic mass divergence at large electron-phonon cou- 
pling strength |44| . which is superficially very similar to 
the quasiparticle effective mass divergence we find in the 
current interacting electron problem. A key difference be- 
tween the polaronic problem and current interacting elec- 
tron problem is the fact that the polaronic case deals with 
a single electron strongly coupled to a bosonic bath (i.e. 
phonons) whereas the current problem involves infinite 
fermionic degrees of freedom (i.e. this is a many-electron 
problem). Thus, the polaronic effective mass divergence 
(the so-called self-trapping of the polaron) , being a one 



electron problem, cannot by definition involve any phase 
transition whereas the current problem being a problem 
with infinite degrees of freedom does (at least) allow the 
possibility of a phase transition. Whether the quasipar- 
ticle effective mass divergence reported in this work is 
indeed a phase transition or not remains to established 
through further work. 

This work is supported by the ONR, the LPS, and the 
NSF. 

Six months after the submission of our 
manuscript for publication and its posting on 
the archive (cond-mat/0312565), a recent preprint 
(cond- mat/0406676) has appeared which verifies our 
finding of the quasiparticle effective mass divergence 
in low density interacting electron liquids at a critical 
density. In this preprint (cond-mat/0406676), the quasi- 
particle effective mass divergence is found to occur even 
within many-body local field approximations which are 
more sophisticated than the simple RPA and Hubbard 
approximation many-body theories utilized in our work. 
The theoretical finding of a quasipasrticle effective 
mass divergence in interacting electron liquids by two 
different groups employing different approximation 
schemes give us confidence that the mass divergence is 
a true many-body phenomenon, and not an artifact of a 
specific approximation scheme. 
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